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Abstract 

A compact approximation formula for the mass of the lightest neutral CP-even 
Higgs boson, m^, in the Minimal Supersymmetric Standard Model (MSSM) is derived 
from the diagrammatic two-loop result for up to (D(aa s ). By analytically expanding 
the diagrammatic result the leading logarithmic and non-logarithmic as well as the 
dominant subleading contributions are obtained. The approximation formula is valid 
for general mixing in the scalar top sector and arbitrary choices of the parameters in 
the Higgs sector of the model. Its quality is analyzed by comparing it with the full 
diagrammatic result. We find agreement with the full result better than 2 GeV for 
most parts of the MSSM parameter space. 
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1 Introduction 



One of the most striking phenomenological implications of Supersymmetry (SUSY) is the 
prediction of a relatively light Higgs boson, which is common to all supersymmetric models 
whose couplings remain in the perturbative regime up to a very high energy scale [0]. The 
search for the lightest Higgs boson thus allows a crucial test of SUSY, and is one of the 
main goals at the present and the next generation of colliders. A precise knowledge of the 
dependence of the mass rrih of the lightest Higgs boson on the relevant SUSY parameters is 
necessary for a detailed analysis of SUSY phenomenology at LEP2, the upgraded Tevatron, 
and also for the LHC and a future linear e + e~ collider, where a high-precision measurement 
of rrih might become possible. 

In the Minimal Supersymmetric Standard Model (MSSM) 0, at the tree level the mass 
of the lightest Higgs boson is restricted to be smaller than the Z-boson mass Mz- This 
bound, however, is strongly affected by radiative corrections, resulting in an upper bound 
of about 135 GeV PHIS]]. Results beyond one-loop order have been obtained using several 
different approaches: a Feynman-diagrammatic calculation of the leading QCD corrections 
has been performed f§-10|; renormalization group (RG) methods have been applied in order 
to obtain leading logarithmic higher-order contributions |TTHl3[; the leading two-loop QCD 
corrections have been calculated in the effective potential method [jl4l, [l5f|. Until recently 
phenomenological analyses have been based either on RG improved one-loop calculations |p"T| - 
|l3l or on the complete Feynman-diagrammatic one- loop on-shell result |5|-|7|. Their numerical 
results, however, differ by up to 20 GeV. Recently the Feynman-diagrammatic result for the 
dominant contributions in 0(aa s ) to the masses of the neutral CP-even Higgs bosons has 
become available ||. By combining these contributions with the complete one-loop on-shell 
result ||, the presently most precise result for rrih based on diagrammatic calculations is 
obtained |||n|. In comparison with the results obtained by RG methods good agreement 
is found in the case of vanishing mixing in the scalar top sector, while sizeable deviations 
which can exceed 5 GeV occur when mixing in the t-sector is taken into account P,[l0[|. 

The Feynman-diagrammatic two-loop result for rrih, however, is very lengthy, making 
the evaluation of the Higgs-boson masses in this approach relatively slow. This could limit 
the applicability of this result e.g. in Monte Carlo simulations. In the present paper we 
derive, by means of a Taylor expansion, a short analytical approximation formula from 
the diagrammatic two-loop result up to 0(aa s ) The purpose of this is not only to 

provide a compact analytical expression for rrih suitable for a very fast numerical evaluation 
without losing too much of accuracy, but also to isolate the most important contributions, 
thus allowing a better qualitative understanding of the source of the dominant corrections. 
The compact approximation formula contains, besides the relevant parts of the one-loop 



contributions, partly taken over from Ref. [O], the leading logarithmic and non- logarithmic 



two-loop corrections for general mixing in the t-sector, the leading Yukawa corrections fll^p: 
and leading QCD contributions beyond 0(aa s ). 

The values for rrih obtained from the approximation formula are compared with the full 
result ||[n|. The dependence on the various MSSM parameters from the stop sector, the 
Higgs sector and the chargino-neutralino sector is analyzed. We find that the approximation 
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formula agrees with the full result better than 2 GeV for most parts of the MSSM parameter 
space. 

The paper is organized as follows: In Sect. 2 a Taylor expansion of the diagrammatic 
two-loop result up to 0(aa s ) is performed and a compact approximation formula is derived. 
The accuracy of the approximation formula is discussed in Sect. 3, where also the location of 
the maximal values of rrih, depending on the mixing in the t-sector, is analyzed analytically. 
In Sect. 4 we give our conclusions. 



2 The compact analytical formula 

2.1 The scalar top sector of the MSSM 

In order to fix our notations and to explain the approximations employed in this sector, we 
first list the conventions for the MSSM scalar top sector: the mass matrix in the basis of the 
current eigenstates tjj and is given by 

KA 2_ (Ml+mt + cos 2/?(±-§4 / )Mf m t M t LR \ 

7W * ~ ^ m t M t LR M? r + m\ + § cos 2(3 s 2 w M 2 z ) ' W 

where 

m t M t LR = mt{A t -n cot (3 ) . (2) 
Neglecting the numerically small contributions proportional to M§ and setting^ 

M- tL = M~ tR := mg, M 2 := m\ + m 2 (3) 

one arives at 

. 2 _ / M| m t M t LR \ 

Mi ~ \ m t M t LR M 2 S ) ■ [A) 

Diagonalizing the t-mass matrix (^) yields the mass eigenvalues mf , and the t-mixing 
angle which relates the current eigenstates to the mass eigenstates: 



?2 J \ — sin 9i cos Qi J I tn 



(5) 



In the above approximation, which we will use throughout the rest of the paper, the 
t-masses and the mixing angle are given by 

ml = M 2 -\m t M t LR \=M 2 (l-At), 

ml = M 2 S + \m t M t LR \= M 2 (1 + A f ) , (6) 



f for M t LR < 



Ot = { for M t LR = , (7) 
for M, LiJ > 



1 Later we will also discuss the case Mi L ^ ^t R i ^ or which the same formalism applies, but with a 
different definition for Ms, see eq. (|l9|) . 



with 

* Ml ml + ml ' K ' 

From the above definition it follows that < A t ~ < 1, otherwise the t-mass matrix eq. (§) 
would have a negative eigenvalue. 



2.2 Calculation of the mass of the lightest Higgs boson 

Here we only give a very brief outline of the calculation of the mass of the lightest neutral CP- 
even Higgs boson in the MSSM; for notations and a detailed description see Refs. pHlD|. We 
focus on the different steps of approximations made in order to derive a compact analytical 
expression from the full diagrammatic result of 0(aa s ). 

At the tree level, the mass matrix of the neutral CP-even Higgs bosons in the basis of 
the weak eigenstates 0i, 2 can be expressed in terms of M z , Ma (the mass of the CP-odd 
Higgs boson) and tan (3 = v 2 /v i (the ratio of the vacuum expectation values of the two Higgs 
doublets) as follows: 



M, 



2,tree _ / Kl^fa 




Higgs - \ m 2 m 2 



<j>2 



Ml sin 2 (3 + M| cos 2 (3 -{M\ + M§) sin (3 cos (3 
-(Ml + Ml) sin (3 cos (3 M 2 A cos 2 (3 + Mf sin 2 (3 



(9) 



Taking into account higher-order corrections, the Higgs-boson masses can to a good 
approximation be obtained by diagonalizing the matrix 

HlggS [ml^-E^M m 2 2 -S, 2 (0) ) ' 1 } 

where the S s (0) (s = <f)i,4'i4>2, < f ) 2) denote the renormalized Higgs-boson self-energies (in 
the 0i,02 basis). As a first step of approximation, the momentum dependence, which is 
numerically rather small, has been neglected in the E s (p 2 ). 

The mass of the lightest Higgs boson receives contributions from all sectors of the MSSM, 
but not all are numerically of equal relevance. In order to derive a compact analytical 
expression, we make the following further approximations: 

The contributions from the t, t-sector up to the two-loop level: The contributions 
arising from the t, t-sector can be written as 

S,(0) = Si 1 )(0)+Si 2) (0), 5 = 01,0102,02- (11) 

Here the E^(0) denote the one-loop contributions of the t, t-sector to the renormalized 
Higgs-boson self-energies. Their explicit form (including also the momentum depen- 
dence) can be found e.g. in Ref. ||. The S^O) denote the two-loop contributions 
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from the t, t-sector at zero external momentum from the Yukawa part of the theory 
(neglecting the gauge couplings) 0. We first consider the dependence of the mass 
shift Am\ (ml = m^ tme + Am^) on tan/? and fi. For the leading contributions to 
Ami from the Yukawa part of the theory the dependence on tan/5 and \i drops out 
in the limit Ma 3> Mz-f\ This holds both for the contributions in O(o) and 0(aa s ). 
Since tan/3 and [i enter only via non- leading corrections in Am\, the dependence on 
them is relatively mild. We therefore use the approximation fi = in the £ s (0). We 
furthermore extract a common prefactor (1/ sin 2 /5) and set otherwise sin/5 = 1 in the 
non-logarithmic one-loop contributions, while the full dependence on sin f3 is kept in 
the logarithmic one-loop and the two-loop contributions. Since the variation of m^ 
with m g is ±2 GeV at most, see Ref. 0], we have eliminated the dependence of the 
S^O) on the gluino mass by setting 



m g = M SUSY = yjM% - ml (12) 

where in the case M^ L = Mi R = m g the SUSY scale is given by Msusy = rn-q- 

As the main step of our approximations, we have performed a Taylor expansion in 
of the S s (0) by inserting eq. (Q) for the t-masses and eq. (0) for the t-mixing angle. For 
the one-loop correction we have expanded up to 0(A~); all three renormalized one- 
loop Higgs-boson self-energies give a contribution. In the one-loop self-energies we have 
kept terms up to 0(M^/mf) while terms of (9(M|/M|) have been neglected. We have 
checked that the numerical effect of the latter terms is insignificant. Concerning the 
two-loop self-energies, the expansion has been carried out up to O(A-). With the above 

described approximations only E?^ gives a non-zero contribution. The Taylor expanded 
self-energies have then been inserted into the Higgs-boson mass matrix eq. ([101) . 

The one-loop contributions from the other sectors: For the one-loop corrections from 
the other sectors of the MSSM we use the logarithmic approximation given in Ref. ||. 
In this approximation the scale of the soft SUSY-breaking parameter in the gaugino 
sector, M, is chosen as M = M SUS y, where M S usy is defined as in eq. ([TJ]). Besides 
Msusy, this contribution is parameterized in terms of Ma, the mixing angle /3, and 
the SM parameters, see eq. (^). A higher accuracy of the non-leading one- loop con- 
tributions in the approximation formula can be achieved by including further terms of 
the one-loop logarithmic approximation given in Ref. |13| . 



Corrections beyond 0(aa s ): Leading contributions beyond 0(aa s ) have been taken into 
account by incorporating the leading two-loop Yukawa correction of 0(G 2 F m^) |T2|,|T6|j 
and by expressing the t, t-contributions through the MS top-quark mass 

m t = m t (m t ) w — ™\ — - (13) 
1 + j^a s (m t ) 



instead of the pole mass m t . This leads to an additional contribution in 0(aa 



2 Note that M t Li? is treated as a free parameter in eq. (|^) and therefore does not depend on /i. 
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2.3 The analytical approximation formula for rrih 
2.3.1 The case for general Ma 



With the approximations described above we obtain the following contributions to the renor- 
malized Higgs-boson self-energies from the t, t-sector (expressed in terms of the top-quark 
pole mass) at one-loop order: 



4>) 



G F y/2 



71" 



Mf A cos 2 ft log 



m t 



41(0) 



M 2 Z cot ft 



7T 

G F y/2 



■-m 2 t + Mf A sin 2 ft 



M 



it 2 8 sin 2 /? 



m 



z 11 M| 
10 m£ 



+ 
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+ 



+ 



12-6 



Mf o „ Mi 4 „ 
z sin 2 /? +8^Asin 4 /? 



m 2 



4 

m? 



Hi, 



(14) 
(15) 



with 



{M t LR j 
Ml 

(M t LR 

~mT 

Mf 



A 



2 

-12 + 4— + 6— 
' M|j 



M 2 



m 2 



+ 



(M^y 

M| 



Mf 



3 m 2 
5M| 



3 m\ 
7 Ml 



12 mf m. 
L + 9 - 

5 Ml Ml 



12 m\ 

Tm| 



+ 



3 

2MJ 



s 



w 



+ 



9 



1 - 



Mi 



14' 



Mf 



(16) 



(17) 



We have verified that the logarithmic terms in eqs. (|14|)-(|16|) agree with the ones given in 
Ref. §. 

The two-loop contributions read: 



4>) 
(0) 



,(2) 



sg(o) 



o , 
o , 

G F \/2 a s 



m. 



7T 



Ti sin ft 



m; 



mi 



31og z hr^ -6 log ~k 



Ml 



Ml 



6 



M t LR 



Mc 



(Mt R )' 
Ml 



mi 



log 



Ms has to be chosen according to 
^jm\ + m 2 t 



Mc 



Mr 



3 (M t LR X 
+ 4" 



Mj 



Mr 
tR 



Mf L Mf R + m 2 t (Ml + MfJ + m 



5 



M lh ± Mr 



(19) 



tR 



The last formula requires some explanation: when performing the expansion in Aj, it was 
assumed that Mi = M~ tR . Thus our result contains only one soft SUSY-breaking scale Mg. 
For the case Mi L ^ Mf R , M s is chosen to reproduce the argument of the leading log correctly. 
According to eq. ( |i~9l) the log yields in both cases 



Mth) = lo g(zr^-) +0(Ai) (20) 



M SJ Wl^fe, 



in agreement with eq. 

The one-loop contributions from the other sectors of the MSSM are not listed here but 
can be found in Refs. P, 13|. The combined self-energies have to be inserted into eq. (|T0D. 



In order to incorporate the leading QCD corrections beyond 0(aa s ), the top-quark pole 
mass should be replaced by the MS top-quark mass eq. fll3|) in the two-loop contribution, as 
described above. Diagonalization of eq. ( |10|) yields the square of the masses of the neutral 



CP-even Higgs bosons. In order to incorporate leading contributions in 0(a 2 ), the leading 
Yukawa correction, see eq. (p8|) below, can be added. 



2.3.2 The case M A > M z 

The diagonalization of the mass matrix eq. (|I~0| ) in the evaluation of m\ incorporates con- 
tributions that are formally of higher order but are non-negligible in general. For large Ma 
these higher-order contributions are suppressed by inverse powers of Ma- Therefore it is 
possible in this case to perform an expansion in the loop order, leading to a very compact 
formula for ml of the form 



ml 



2, tree . \ 2,a,t/t , « 2, a, rest . \ 2,aa, 



+ Am h . 



(21) 



At the two-loop level, only the <p2 self-energy contributes in our approximation, yielding 
the term 



with E?^(0) from eq. 



Am^ = -sin a 0Eg(O) (22) 

In contrast to eqs. (fn])-(|TE|) we give here an expression for rrih 
in which the MS top-quark mass is used everywhere. Inserting eq. ([13]) into E^J(O) yields 
additional contributions to Am^' Qtts , see eq. (BTj) . From eq. ( ^2] ) it becomes also clear that 
the 13 dependence of the leading two-loop contribution drops out, as mentioned above. 

The tree-level and one-loop contributions of eq. ([H]) are given by 



m 



2, tree 
h 



Am 2 h a > t/E 



M\ + M 2 Z - J (Ml + M|) 2 - AM 2 Z M\ cos 2 2(5 



(23) 



G F V2 



m. 



71 z 



log ( S ) { A ~ V^r cos 2/3 - cos 2 2(3 



' 9 

1717 



Ml 



4 mf 



M 2 ( 3 M 2 \l 1 

z cos 2 p cos 2(3 (6 + -^f (1 -4 sin 2 /?) - ^8A cos 2(3 sin 2 /? 



M\ 



2 m 2 



6 



+ 



+ 



mf 



\\M% (M t LR )' 
80 mf + 



Ml 



(m^) 4 / 


1 1 mf 


3 m 4 


Ml [ 


8 ' 2M| 


8M| 



3 

) + 



lMf 
2 mf 

(M t LR 
~Mj 



3mf\ 
4Ml) 



3 mf 
' 40 MJ 



3 m 4 
10 M| 



4MfJ 



+ 



A 2, o, rest 



Ml 
G F V2Ml 
~V 2 2A 



3 mf 
56Mf 



+ 



3 mf 
14 Mf 



3 7nf\ 
16MfJ 



1 + 4 



Mf 



cos 2 (3 cos 2/9 



(24) 



log 



. Ml 



- 



12N r 



Ml 



ni 



6N C cos 2/3 — | + cos 2/? (P 6 + P f ) 



Mf 



+ (P g + P2h) (sin 4 (3 + cos 4 /? ) - 2 sin 2 cos 2 (3 (PL + P' 2H ) + 



QN c mj 
Ml 



{Mt R ) 2 
M 2 

IVI SUSY 



( 



M b LR \ 



12 Mg USY 



Ml 

— ^2 cos 2(3 
2mf 



M b LR 



M 2 

1V1 SUSY 



+ 



A\ - /i 2 tan 2 (3 
3 Mg USY 



-log 



(cos 4 (3 + sin 4 /3)P 2H -2 cos 2 /3 sin 2 /3 P 2 ' H - P 1H 



where A is defined as in eq. (|T7D, and 



M 



p, 



Ah — /i tan /3 

^(1 + 43,4, + 8g 2 s ^; 
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PlH 
P2H 



bJ S$y)+N g (2-4s$y + 8s$y) 
2 o«„4 



1 

W 



N c (N g -l)(2-is 2 w + S(Q 2 t 

-44 + 106^-62^, P; = 10 + 344,-26^ 
-9 cos 4 2(3 + (1 - 2^ + 25^) cos 2 2(3 

-10 + 24 -24, p^ = 8-224,- 

Qt = \, G& = -~, N c = 3, N g = 3. 
The dominant two-loop contribution of 0(aa s ) to ml reads 



Am 



2,aa„ 
7l 



m. 



7T Z 



7T 



/ ml 



+ 2 log 



+ 



Mf k 

n(M t LR )' 



3 log 



12 



Ml 



1 + 4 



Mf 
Mf 



cos 2 (3 cos 2/? 



,M f LK 



M, 



(25) 



(26) 



(27) 



For the leading two-loop Yukawa correction we use the result given in [^2, 16fl: 

9 



Am 



2,a z 



16vr 4 F 1 



xt + r 



(28) 



2 2 

mf — mi 

i to ti • 2 

H — L sm z 



2m^ 



2% log 



1 / m- m~ 

2 log Hsr 1 



(29) 



Here M5 is chosen in analogy with eq. ([19]), where the top-quark pole mass m t has to be 
replaced by the MS top-quark mass m t . 

In the contributions of the t,t-sector at one- loop and two- loop order, eqs. ( P^[ ) an d (|2"TD, 
we have included correction factors of 0(M^/M\). In this way the compact formula eq. (|2l|) 
gives a reliable approximation for Ma values down to at least Ma = 200 GeV. 

The contribution of 0(aa s ) given in eq. ( p7|) can be compared with analytical formulas 
derived via the two-loop effective potential approach for the case of no mixing in the t- 
and via RG methods 0,0. The leading term ~ log 2 (m 2 /Mj) agrees with 
The subleading term for vanishing t-mixing ~ log(m 2 /Mj) 
agrees with the result of the two-loop effective potential approach [|T3[ and the result of the 



sector im 



the results in Refs. 12-14 



two-loop RG calculation |L3|,|14||, but differs from the RG improved one-loop result [|l"2||T3| 
The term ~ log(m 2 / Mj)(M t LR / Ms) 2 for non- vanishing t-mixing differs from the result given 
in Ref. [P2], 13fl. All other terms of 0(aa s ) are new. The term ~ Mf R /Ms shows that the 
result for mn is not symmetric in ±M^ R . The good numerical agreement with the RG results 
in the case of no mixing in the t-sector can qualitatively be understood by noting that in 
the no-mixing case the leading term in both approaches agrees, while for the corrections 
proportional to powers of M^ R /Ms deviations occur already in the leading contribution. 

When comparing the results of the diagrammatic on-shell calculation with the RG results 
in terms of the parameters Ms and M^ R it should be noted that, due to the different 
renormalization schemes employed, the meaning of these (non-observable) parameters is not 
precisely the same in the two approaches starting from two-loop order (see the discussion in 
Ref. U). A more detailed comparison of our results with those obtained via RG methods 
will be performed in a forthcoming publication. 



2.4 Implementation into FeynHiggs 

The formulas given in the previous section have been implemented into the Fortran code 
FeynHiggs [O, thus allowing a direct comparison between the full result described in 



Refs. |||T0| and the approximation formula (^T|) . 

We also provide the Fortran code FeynHiggsFast in which only the formula for the com- 
pact approximation of m^ is implemented, thus allowing an approximate but much faster 
evaluation of the Higgs-boson mass m^. This program is shorter by a factor of 50 and faster 
by a factor of 3 x 10 4 with respect to FeynHiggs] it needs only 2 x 10~ 5 seconds on a Sigma 
station (Alpha CPU, 600 MHz) for the evaluation of mh for one set of parameters. The qual- 
ity of the prediction of m^ is reasonably good for a large part of the MSSM parameter space, 
as will be discussed in the next section. Into FeynHiggsFast we have also implemented 
the calculation of the MSSM contributions to Ap |T3. Here the corrections arising from 



S 



t/b- loops up to 0(aa s ) have been taken into account, neglegting only the gluino-exchange 
contribution which is very lengthy and vanishes for large m g . The p-parameter can be used 
as an additional constraint (besides the experimental bounds) on the squark masses. A 
value of Ap outside the experimentally preferred region of Ap SUSY < 1 x 10~ 3 [|19| indicates 
experimentally disfavored t- and 6-masses. 

Both Fortran codes can be obtained via the WWW page 
trttp : / /www-itp . physik . uni-karlsruhe . de/f eynhiggs . 



3 Discussion of the compact expression for rrih 

In this section we analyze the quality of our compact approximation formula for rrih with 
respect to the full calculation, which contains the full diagrammatic one-loop contribution || , 
the complete leading two-loop corrections in 0(aa s ) pHIOll, and the two contributions beyond 
0(aa s ), see eq. ( |13|) and eq. (|28|) . We use here the approximation formula for general Ma, 
as discussed in Sec. |2.3.1| . As mentioned above, in the region Ma ^ 200 GeV the compact 



formula given in Sec. |2.3.2| yields equally good results. 



For tan f3 we restrict ourselves to two typical values which are favored by SUSY-GUT 
pOp : tan/3 = 1.6 for the SU(5) scenario and tan/3 = 40 for the SO (10) sce- 



scenarios 



nario. Other parameters are M z = 91.187 GeV, M w = 80.39 GeV, G F = 1.16639 x 
10~ 5 GeV -2 , a s (m t ) = 0.1095, m t = 175 GeV, and mj = 4.5 GeV. In the numerical eval- 
uation we have furthermore chosen A t = A^. The parameter M appearing in the plots is 
the SU(2) gaugino mass parameter, it enters in the full result only, see the discussion of 



the one-loop contributions in Sect. |2.2| . p is the Higgs- mixing parameter. If not indicated 
differently, we have chosen M = rriq and fi = —rriq. 

In Fig. [I] we show rrih as a function of M^ R /rriq for rriq = 200,500,1000 GeV and 
Ma = 500 GeV. A maximum for rrih, evaluated with the full formula as well as with the 
approximation formula, is reached for about M^ R /rriq « ±1.9 in the tan/? = 1.6 scenario 
and in the tan/? = 40 scenario. This case we refer to as 'maximal mixing'. A minimum is 
reached around Mf R /rriq ~ which we refer to as 'no mixing'. In general the approximation 
differs from the full result by less than 2 GeV up to \M^ R /rriq\ < 2. For larger \M^ R \ size- 
able deviations occur, which become very large for \M^ R /rriq\ > 2.5. The reason is that the 
expansion parameter becomes rather large and approaches 1 in this region. The effect 
of the non-logarithmic terms in the two-loop contribution in 0(aa s ) reaches up to about 
5 GeV for maximal mixing. 

The location of the extrema can be understood analytically as follows: taking into account 
only the leading one-loop corrections from the t,i-sector for the calculation of rrih, one can 
easily compute the values of M^ R / rriq for which rrih reaches a maximum or a minimum. The 
well known results are 

i \/6 (maximum) 

M LR 

= I (minimum) . (30) 

\ — v 6 (maximum) 
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Figure 1: as a function of M^ R /rrig, calculated from the full formula and from the 

approximation formula for Ma = 500 GeV, m- g = 500 GeV and tan/3 = 1.6 or 40. 



2,aa a 



Taking into account the new two-loop corrections Am ft ' ° , given m eq. 
of the extrema in eq. (|30|) receive a shift, yielding up to 0(aa s 



7f), the positions 



M t LR 
ma 



v 7T 
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1 + 3v/6- v^log 



Mi 



+1.92 for M s = 1000 GeV: max.) 



—0.07; minimum) 
-1.85 for M s 



1000 GeV; max.) 

(31) 

The maxima are shifted to smaller absolute values of M^ R /rrig, the minimum is shifted to a 



Mi 



slightly negative value (see the discussion in Refs. 

Fig. |2| shows the dependence of on rrtq for the cases of no mixing and maximal mixing, 
and we have set Ma = 500 GeV. Very good agreement is found in the no-mixing scenario 
as well as in the maximal-mixing scenario, the deviation lies below 2 GeV. 

The dependence on Ma is shown in Fig. RL The quality of the approximation is typically 
better than 1 GeV for the no-mixing case and better than 2 GeV for the maximal-mixing 
case. Only for very small (and experimentally already excluded) values of Ma a deviation 
of 5 GeV occurs. The peaks in the plot for tan/3 = 1.6 in the full result are due to the 
threshold Ma = 2 rrit in the one-loop contribution, originating from the top-loop diagram in 
the A self-energy. This peak does not occur in the approximation formula and can thus lead 
to a larger deviation around the threshold. 



In deriving the leading and subleading corrections we have set m- Q 



Ml 



m 



t ■ 



thus 



eliminating this additional scale. A variation of rrig in the full result directly corresponds to 
a shift relative to the approximation formula. As it was shown in Ref. |10|| , this deviation is 
negligible for the no-mixing case and lies within ±2 GeV in the maximal-mixing case. It was 
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Figure 2: as a function of m,, calculated from the full formula and from the approxi- 
mation formula for Ma = 500 GeV, m~ g = 500 GeV and tan (3 = 1.6 or 40. 
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Figure 3: rrih as a function of Ma, calculated from the full formula and from the approxi- 
mation formula for irtq = 1000 GeV, rrig = 500 GeV and tan/? = 1.6 or 40. 



also analyzed in Ref. [Mj that the variation of with the SU{2) gaugino mass parameter 
M(= M 2 ) and with the Higgs-mixing parameter [i is relatively weak. We have found only 
small deviations of the approximation formula with respect to the full result, not more than 
2 GeV for small values of M, and 1.5 GeV for small values of \x. 

Finally we consider the case when M$ L 7^ Mi R , which is shown in Fig. [|. Here also the 
case Mt R = Mr is depicted, which we refer to as 'moderate mixing'. 'No mixing' here 
corresponds to M^ R = 0, 'maximal mixing' corresponds to the choice M^ R = 2 M^ L . For 
fixed Mi R this leads to a t-mass below the experimental lower bound for nearly the whole 
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Figure 4: as a function of Mf R or Mf L , calculated from the full formula and from the 
approximation formula for varied M {r (M^), M h (M- tR ) = 300 GeV, M A = 500 GeV,m § = 
500 GeV and tan (3 = 1.6 or 40. 



parameter space; thus we have omitted this case in this scenario (bottom of Fig. Here 
we have made use of eq. flT§p. Agreement better than 2 GeV is found for all scenarios. 

4 Conclusions 

By means of a Taylor expansion of the diagrammatic two-loop result we have derived the 
leading logarithmic and non-logarithmic as well as subleading contributions to the lightest 
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Higgs-boson mass m h up to 0(aa s ). This result has been incorporated into a compact 
analytical formula for m\, which can easily be implemented for numerical evaluation. It 
contains the tree-level expression, leading and subleading one-loop contributions, the leading 
logarithmic and non-logarithmic two-loop contributions, and two further corrections beyond 
0(aa s ). This formula is valid for general mixing in the scalar top sector and arbitrary choices 
of the parameters in the Higgs sector of the model. It has been included into the Fortran codes 
FeynHiggs and FeynHiggsFast, thus allowing a very fast and reasonably accurate evaluation 
of mh- The codes are available via the WWW page 
|http : //www-itp . physik . uni-karlsruhe . de/f eynhiggs . 

The approximation formula has been compared with the full diagrammatic result, in 
which the complete one-loop and the dominant two- loop corrections of 0(aa s ) are taken 
into account without Taylor expansion. Good agreement, better than 2 GeV, is found for 
most parts of the MSSM parameter space. Larger deviations, exceeding a few GeV, can occur 
for large mixing in the t-sector, |M t / 'rrtq\ > 2. The effect of the corrections of 0(aa s ) of 
shifting the maxima for rrih towards smaller values of \M^ R /rriq\ can easily be read off from 
the compact formula by analytically determining the extrema of rrih- 
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